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On the 2:1 Orbital Resonance in the HD 82943 Planetary System^ 

Man Hoi Lee^, R. Paul Butler^, Debra A. Fischer‘S, Geoffrey W. Marcy^, and Steven S. Vogt® 

ABSTRACT 

We present an analysis of the HD 82943 planetary system based on a radial velocity 
data set that combines new measurements obtained with the Keck telescope and the 
CORALIE measurements published in graphical form. We examine simultaneously the 
goodness of fit and the dynamical properties of the best-fit double-Keplerian model as 
a function of the poorly constrained eccentricity and argument of periapse of the outer 
planet’s orbit. The fit with the minimum is dynamically unstable if the orbits are 
assumed to be coplanar. However, the minimum is relatively shallow, and there is a 
wide range of fits outside the minimum with reasonable Xu- assumed coplanar 

inclination i = 30° (sin i = 0.5), only good fits with both of the lowest order, eccentricity- 
type mean-motion resonance variables at the 2:1 commensurability, 6 i and 62 , librating 
about 0° are stable. For sini = 1, there are also some good hts with only 9i (involving 
the inner planet’s periapse longitude) librating that are stable for at least 10® years. 
The libration semiamplitudes are about 6° for 61 and 10° for 62 for the stable good 
fit with the smallest libration amplitudes of both 81 and 62 - We do not find any good 
fits that are non-resonant and stable. Thus the two planets in the HD 82943 system 
are almost certainly in 2:1 mean-motion resonance, with at least 9i librating, and the 
observations may even be consistent with small-amplitude librations of both 9i and 02 - 


1. INTRODUCTION 

The hrst pair of extrasolar planets suspected to be in mean-motion resonance was discovered 
around the star GJ 876, with the orbital periods nearly in the ratio 2:1 (Marcy et al. 2001). A 
dynamical fit to the radial velocity data of GJ 876 that accounts for the mutual gravitational 
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interaction of the planets is essential because of the short orbital periods 30 and 60 days) 
and large planetary masses [minimum combined planetary mass relative to the stellar mass (mi + 
'ni 2 )/mo ss 0.0074], It is now well established that this pair of planets is deep in 2:1 orbital 
resonance, with both of the lowest order, eccentricity-type mean-motion resonance variables, 

01 = Ai — 2 A 2 -|- tJ7l (1) 


and 

02 = Ai — 2 A 2 -|- W 2 -, (2) 

librating about 0° with small amplitudes (Laughlin & Chambers 2001; Rivera & Lissauer 2001; 
Laughlin et al. 2005). Here Ai ^2 are the mean longitudes of the inner and outer planets, respectively, 
and wi ^2 are the longitudes of periapse. The simultaneous librations of 0i and 02 about 0° mean 
that the secular apsidal resonance variable, 

0SAR = Wi - ^2 = 9l — 02, (3) 


also librates about 0° and that the periapses are on average aligned. The basic dynamical properties 
of this resonant pair of planets are not affected by the recent discovery of a third, low-mass planet 
on a 1.9-day orbit in the GJ 876 system (Rivera et al. 2005). 

R is important to confirm other suspected resonant planetary systems, as the GJ 876 system 
has shown that resonant planetary systems are interesting in terms of both their dynamics and their 
constraints on processes during planet formation. The geometry of the 2:1 resonances in the GJ 
876 system is different from that of the 2:1 resonances between the Jovian satellites lo and Europa, 
where 0i librates about 0° but 02 and 0sar librate about 180°. For small orbital eccentricities, the 
lo-Europa configuration is the only stable 2:1 resonance configuration with both 0i and 02 librating. 
For moderate to large eccentricities, the lo-Europa configuration is not stable, but there is a wide 
variety of other stable 2:1 resonance configurations, including the GJ 876 configuration (Lee & Peale 
2002, 2003a; Beauge et al. 2003; Ferraz-Mello et al. 2003; Lee 2004). The resonances in the GJ 
876 system were most likely established by converging differential migration of the planets leading 
to capture into resonances, with the migration due to interaction with the protoplanetary disk. 
While it is easy to establish the observed resonance geometry of the GJ 876 system by convergent 
migration, the observational upper limits on the eccentricities (ei < 0.31 and 62 < 0.05) require 
either significant eccentricity damping from planet-disk interaction or resonance capture occurring 
just before disk dispersal, because continued migration after resonance capture could lead to rapid 
growth of the eccentricities (Lee &: Peale 2002; Kley et al. 2005). Hydrodynamic simulations of the 
assembly of the GJ 876 resonances performed to date do not show significant eccentricity damping 
from planet-disk interaction and produce eccentricities that exceed the observational upper limits, 
unless the disk is dispersed shortly after resonance capture (Papaloizou 2003; Kley et al. 2004, 
2005). 

HD 82943 was the second star discovered to host a pair of planets with orbital periods nearly 
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in the ratio 2:1. The discovery of the first planet was announced in ESO Press Release 13/00^ and 
the discovery of the second, inner planet was announced in ESO Press Release 07/01^. Mayor et 
ah (2004) have recently published the radial velocity measurements obtained with the CORALIE 
spectrograph on the 1.2-m Euler Swiss telescope at the ESO La Silla Observatory in graphical form 
only. Unlike GJ 876, a double-Keplerian fit is likely adequate for HD 82943, because the orbital 
periods are much longer (~ 220 and 440 days) and the planetary masses are smaller [minimum 
(mi + m 2 )/mo K. 0.003], and the mutual gravitational interaction of the planets is not expected 
to affect the radial velocity significantly over the few-year time span of the available observations. 
Mayor et al. (2004) have found a best-fit double-Keplerian solution, and its orbital parameters are 
reproduced in Table 1. 

Eerraz-Mello et al. (2005) have reported simulations of the Mayor et al. (2004) best-fit solution 
that are unstable, but they assumed that the orbital parameters are in astrocentric coordinates 
and they did not state the ranges of orbital inclinations and starting epochs examined. The fact 
that the orbital parameters obtained from multiple-Keplerian fits should be interpreted as in Jacobi 
coordinates was hrst pointed out by Lissauer & Rivera (2001) and was derived and demonstrated 
by Lee & Peale (2003b). We have performed direct numerical orbit integrations of the Mayor et 
al. (2004) best-fit solution, assuming that the orbital parameters are in Jacobi coordinates and 
that the orbits are coplanar with the same inclination i from the plane of the sky. When we 
assumed that the orbital parameters correspond to the osculating parameters at the epoch JD 
2451185.1 of the first CORALIE measurement, the system becomes unstable after a time ranging 
from ~ 500 to 4 x 10^ yr for sini = 1, 0.9, ..., 0.5. Since the long-term evolution can be sensitive to 
the epoch that the orbital parameters are assumed to correspond to, we have repeated the direct 
integrations by starting at three other epochs equally spaced between the first (JD 2451185.1) and 
last (JD 2452777.7) CORALIE measurements. All of the integrations become unstable, with the 
vast majority in less than 10® yr. These results and those in Eerraz-Mello et al. (2005, who also 
examined mutually inclined orbits and the effects of the uncertainty in the stellar mass) show that 
the best-fit solution found by Mayor et al. (2004) is unstable. 

It is, however, essential that one examines not only the best-fit solution that minimizes the 
reduced chi-square statistic but also fits with not significantly above the minimum, especially 
if the minimum is shallow and changes slowly with variations in one or more of the parameters. 
Eerraz-Mello et al. (2005) have analyzed the CORALIE data published in Mayor et al. (2004) by 
generating a large number of orbital parameter sets and selecting those that fit the radial velocity 
data with the RMS of the residuals (instead of xV) dose to the minimum value. The stable coplanar 
fits that they have found have arguments of periapse lOi ss 120° and (jJ 2 > 200°, and both 9i and 62 
librate about 0° with large amplitudes.^ On the other hand, Mayor et al. (2004) have noted that 
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it is possible to find an aligned configuration with u)i ~ u )2 that fits the data with nearly the same 
RMS as their best-fit solution, although it is unclear whether such a fit would be stable and/or in 
2:1 orbital resonance. Thus, it has not been firmly established that the pair of planets around HD 
82943 are in 2:1 orbital resonance, even though it is likely that a pair of planets of order Jupiter 
mass so close to the 2:1 mean-motion commensurability would be dynamically unstable unless they 
are in 2:1 orbital resonance. 

In this paper we present an analysis of the HD 82943 planetary system based on a radial 
velocity data set that combines new measurements obtained with the Keck telescope and the 
CORALIE measurements published in graphical form. The stellar characteristics and the radial 
velocity measurements are described in §2. In §3 we present the best-fit double-Keplerian models 
on a grid of the poorly constrained eccentricity, 62, and argument of periapse, of the outer 
planet’s orbit. In §4 we use dynamical stability to narrow the range of reasonable fits and examine 
the dynamical properties of the stable fits. We show that the two planets in the HD 82943 system 
are almost certainly in 2:1 orbital resonance, with at least 9i librating, and may even be consistent 
with small-amplitude librations of both 9i and 02- Our conclusions are summarized and discussed 
in S5. 


2. STELLAR CHARACTERISTICS AND RADIAL VELOCITY 

MEASUREMENTS 

The stellar characteristics of the GO star HD 82943 were summarized in Mayor et al. (2004) 
and Fischer & Valenti (2005). Santos et al. (2000a) and Laws et al. (2003) have determined a stellar 
mass of mo = I.OSMq and I.I8lgQ^M0, respectively. Fischer & Valenti (2005) have measured the 
metallicity of HD 82943 to be [Fe/H] = -1-0.27 which, when combined with evolutionary models, 
gives a stellar mass of I.ISMq. We follow Mayor et al. (2004) and adopt a mass of mo = 1 . 15 M 0 . 

Stellar radial velocities exhibit intrinsic velocity “jitter” caused by acoustic p-modes, turbulent 
convection, star spots, and flows controlled by surface magnetic fields. The level of jitter depends 
on the age and activity level of the star (Saar et al. 1998; Gumming et al. 1999; Santos et al. 
2000b; Wright 2005). From the observed velocity variance of hundreds of stars with no detected 
companions in the California and Carnegie Planet Search Program, Wright (2005) has developed 
an empirical model for predicting the jitter from a star’s B — V color, absolute magnitude, and 
chromospheric emission level. Some of this empirically estimated “jitter” is no doubt actually 
instrumental, stemming from small errors in the data analysis at levels of 1-2 ms“^. Based on the 
model presented in Wright (2005), we estimate the jitter for HD 82943 to be 4.2ms“^, with an 
uncertainty of 50%. 


be measured from any reference direction in the orbital plane. If we choose the ascending node referenced to the 
plane of the sky as the reference direction, zzij are the same as the arguments of periapse coj obtained from the radial 
velocity fit. 
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We began radial velocity measurements for HD 82943 in 2001 with the HIRES echelle spec¬ 
trograph on the Keck I telescope. We have to date obtained 23 velocity measurements spanning 
3.8 yr, and they are listed in Table 2. The median of the internal velocity uncertainties is 3.0 ms“^. 

Mayor et ah (2004) have published 142 radial velocity measurements obtained between 1999 
and 2003 in graphical form only. We have extracted their data from Figure 6 of Mayor et al. (2004) 
with an accuracy of about 0.5 d in the times of observation and an accuracy of about 0.07ms“^ 
in the velocities and their internal uncertainties. The uncertainty of 0.5 d in the extracted time of 
observation can translate into a velocity uncertainty of as much as 2ms“^ in the steepest parts 
of the radial velocity curve (where the radial velocity can change by as much as 4ms“^ per day), 
but the square of even this maximum (2ms“^) is much smaller than the sum of the squares of the 
extracted internal uncertainty (> 4ms“^) and the estimated stellar jitter (4.2ms“^). As another 
check of the quality of this extracted CORALIE data set, we have performed a double-Keplerian 
fit of the extracted data, with the extracted internal uncertainties as the measurement errors in 
the calculation of the reduced chi-square statistic Xu- The resulting best fit has RMS = 6.99 ms“^, 
which is only slightly larger than the RMS (= 6.8 ms“^) of the Mayor et al. best fit, and the orbital 
parameters agree with those found by Mayor et al. and reproduced in Table 1 to better than ±2 
in the last digit. This confirms that the extracted data set is comparable in quality to the actual 
data set. 

In the analysis below, we use the combined data set containing 165 velocity measurements and 
spanning 6.1 yr, and we adopt the quadrature sum of the internal uncertainty and the estimated 
stellar jitter as the measurement error of each data point in the calculation of Xu- 


3. DOUBLE-KEPLERIAN FITS 

As we mention in §1, the mutual gravitational interaction of the planets is not expected to 
affect the radial velocity of HD 82943 significantly over the 6.1-year time span (or about 5 outer 
planet orbits) of the available observations. Thus we fit the radial velocity Vr of the star with the 
reflex motion due to two orbiting planetary companions on independent Keplerian orbits: 

Vr = Ki [cos(u;i + fi) + ei cos cui] -|- K 2 [cos(u;2 + /2) + 62 cos 102 ] + V, (4) 

where V is the line-of-sight velocity of the center of mass of the whole planetary system relative 
to the solar system barycenter, Kj is the amplitude of the velocity induced by the jth planet, 
and Cj, ujj and fj are, respectively, the eccentricity, argument of periapse, and true anomaly of 
the jth planet’s orbit (see, e.g., Lee & Peale 2003b). The zero points of the velocities measured 
by different telescopes are different, and U is a parameter for each telescope (denoted as Vk and 
Vc for the Keck and CORALIE data, respectively). Since the true anomaly fj depends on ej, 
the orbital period Pj, and the mean anomaly Mj at a given epoch (which we choose to be the 
epoch of the first CORALIE measurement, JD 2451185.1), there are five parameters, Pj, Kj, ej, 
ujj, and Mj, for each Keplerian orbit. Thus there are a total of 12 parameters. The fact that 
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the orbital parameters obtained from double-Keplerian fits should be interpreted as Jacobi (and 
not astrocentric) coordinates does not change the fitting function in equation (4) but changes the 
conversion from Kj and Pj to the planetary masses mj and orbital semimajor axes aj: 




2ttG \ rrij sin ij 


m 


I 2/3 


3 




(5) 


Pj = 27r[a|/(Gm')]^/^ ( 6 ) 

where m' = X]i=o ^ Peale 2003b). All of the orbital fits presented here are obtained using 

the Levenberg-Marquardt algorithm (Press et al. 1992) to minimize The Levenberg-Marquardt 
algorithm is effective at converging on a local minimum for given initial values of the parameters. 

We begin by using the parameter values listed in Table 1 as the initial guess and allowing all 
12 parameters to be varied in the fitting. The resulting fit I is listed in Table 3 and compared to 
the radial velocity data in Figure 1. We also show the residuals of the radial velocity data from 
fit I in the top panel of Figure 4. Fit I has x^ = 1-87 for 12 adjustable parameters (or 1.84 if we 
rescale x^ to 10 adjustable parameters for comparison with the 10 parameter fits below) and RMS 
= 7.88 ms“^. However, fit I is dynamically unstable if the orbits are assumed to be coplanar. We 
perform a series of direct numerical orbit integrations of fit I similar to that described in §1 to 
test the stability of the Mayor et al. (2004) best-fit solution (i.e., with sini = 1, 0.9, ..., 0.5 and 
four starting epochs equally spaced between JD 2451185.1 and JD 2452777.7). Figure 2 shows the 
evolution of the semimajor axes, ai and 02 , and eccentricities, ei and 62 , for the integration with 
sini = 1 and starting epoch JD 2451185.1, which becomes unstable after ~ 3200 yr. All of the 
integrations of fit I become unstable in time ranging from ~ 60 to 3200 yr. 

We then explore fits in which one or two of the parameters are fixed at values different from 
those of fit I and the other parameters are varied. These experiments indicate that 62 and U32 are 
the least constrained parameters, as there are hts with 62 and u )2 fixed at values very different from 
those of fit I that have x^ (after taking into account the difference in the number of adjustable 
parameters) and RMS only slightly different from those of fit I. Therefore, we decide to examine 
systematically best-fit double-Keplerian models with fixed values of 62 and 0 J 2 - We start at 62 = 0.22 
and u )2 = 290° (which are close to the values for fit I), using the parameters of fit I as the initial 
guess for the other 10 parameters. Then for each point on a grid of resolution 0.01 in 62 and 5° in 
u) 2 -, we search for the best fit with fixed 62 and u) 2 -, using the best-fit parameters already determined 
for an adjacent grid point as the initial guess. 

Figure 3 shows the contours of x^, RMS, ei, wi, and the planetary mass ratio m\jm 2 for the 
best-fit double-Keplerian models with fixed values of 62 and a; 2 - Fit I at 62 = 0.219 and W 2 = 284° 
is the only minimum in x^j with xt — 1-84 for 10 adjustable parameters. However, the minimum 
is relatively shallow, with Xu ~ ^.0 in a region that includes the entire region with 62 < 0.2 and 
extends to 62 ~ 0.38 near 032 = 130° and 280°. The contours of RMS are similar to those of Xu- 
The RMS minimum (where RMS = 7.87ms“^) at 62 = 0.22 and ijJ 2 = 275° is only slightly offset 



-7 - 


from the Xu minimum, and the region with Xu ~ ^.0 also has RMS < 8.2 ms“^. The contours of 
ei, oji, and mi/m 2 show that these parameters are indeed well constrained, with ei within ±0.05 
of 0.39, u)i within ±15° of 120°, and mi/m 2 within ±0.2 of 1.0 for the fits with x^ < 2.0. 

In Figure 3o we mark the positions of three fits II-IV (as well as fit I) whose dynamical 
properties are discussed in detail in §4. Their orbital parameters are listed in Table 3, and the 
residuals of the radial velocity data from these fits are shown in the lower three panels of Figure 
4. We can infer from the parameters listed in Table 3 that the fits II-IV have nearly identical 
ai(w 0.747AU), a 2 (~ 1.18AU), and m 2 sini(s:i 1.71 Mj, where Mj is the mass of Jupiter) while 
mi sin i ranges from 1.58 Mj for fit II to 1.78 Mj for fit IV. Figure 4 demonstrates visually that 
these hts with Xu only slightly higher than the minimum are comparable in the quality of fit to the 
fit I with the minimum Xu- 


4. DYNAMICAL ANALYSIS 

We begin with some general remarks on the dynamical properties of the best-fit double- 
Keplerian models with fixed 62 and UJ 2 that can be inferred from the contour plots in Figure 
3. There is a wide variety of stable 2:1 resonance configurations with both 9i and 62 librating (see, 
e.g., Lee 2004). For mi/m 2 ~ 1, there are (1) the sequence of resonance configurations that can 
be reached by the differential migration of planets with constant masses and initially coplanar and 
nearly circular orbits, ( 2 ) the sequence with large 62 and asymmetric librations of both 9i and 02 
far from either 0° or 180°, and (3) the sequence with 9i librating about 180° and 02 librating about 
0°. Their loci in the 61-62 plane are shown in Figures 5, 11, and 14, respectively, of Lee (2004). 
The good fits in Figure 3 have 61 ~ 0.39 ± 0.05 and 62 < 0.38, which are far from the loci of the se¬ 
quences 2 and 3. Thus the good fits in Figure 3 cannot be small-libration-amplitude configurations 
in the sequences 2 and 3. Figure 3d shows that the fits with uj 2 ~ 120° also have uji ~ 120°. These 
fits with 0SAR = roi — W 2 ~ 0° currently also have 61 > 0.40 and mi/m 2 < 1.0, and we can infer 
from Figure 5 of Lee (2004) for the sequence 1 that the fits with 62 ~ 0.14 and uj 2 ~ 120° could be 
in a symmetric conhguration with both 0i and 02 (and hence 0 sar) librating about 0° with small 
amplitudes. However, direct numerical orbit integrations are needed to determine whether the fits 
with 62 ~ 0.14 and u )2 ~ 120° are in fact in such a resonance configuration. Finally, Figure 5 of 
Lee (2004) also shows that there are configurations in the sequence 1 with asymmetric librations 
of 01 and 02 in the region 61 0.39 ± 0.05 and 62 < 0.38 if mi/m 2 > 0.95, and we cannot rule out 

the possibility of asymmetric librations for the fits in Figure 3 with mi/m 2 ^ 0.95 without direct 
numerical orbit integrations. 

In order to use dynamical stability to narrow the range of reasonable fits and to determine 
the actual dynamical properties of the stable fits, we perform direct numerical orbit integrations 
of all of the best-fit double-Keplerian models with fixed 62 and uj 2 shown in Figure 3. We use 
the symplectic integrator SyMBA (Duncan et al. 1998) and assume that the orbital parameters 
obtained from the fits are in Jacobi coordinates (Lee & Peale 2003b; see also Lissauer Sz Rivera 



2001). Double-Keplerian fits do not yield the inclinations of the orbits from the plane of the sky. 
We assume that the orbits are coplanar and consider two cases: sini = 1 and 0.5. We perform 
integrations for only one starting epoch, JD 2451185.1, since for each fit we already examine 
neighboring fits with slightly different orbital parameters. The maximum time interval of each 
integration is 5 x 10*^ yr. An integration is stopped and the model is considered to have become 
unstable, if (1) the distance of a planet from the star is less than 0.075 AU or greater than 3.0 AU 
or ( 2 ) there is an encounter between the planets closer than the sum of their physical radii for an 
assumed mean density of lgcm“^. For these planets with orbital periods of about 220 and 440 
days, oi ~ 0.75 AU, 02 ~ 1.2 AU, and their physical radii are about 0.011 and 0.0066 of their Hill 
radii for a mean density of lgcm“^. 

In Figure 5 we present the results of the direct numerical orbit integrations for sini = 1. Figure 
5a shows the contours of the time, tun; at which a model becomes unstable, and Figure 5b shows 
the contours of the semiamplitude A^i. All of the fits that are stable for at least 5 x 10^ yr have 9i 
librating about 0° with A0i < 140°. Figure 5c shows the contours of the semiamplitude A 02 , with 
the thick solid contour indicating the boundary for the libration of 02 about 0°. None of the stable 
fits shows asymmetric libration of 9i or 02 - Since the region with 02 librating is smaller than that 
with 9i librating, there are stable good fits with only 9i librating about 0 °, in addition to stable 
good fits with both 9i and 02 librating about 0°. The fits with 62 ~ 0.14 and a ;2 ~ 130° have small 
amplitude librations of both 9i and 02 about 0°, in good agreement with the estimate of 62 ~ 0.14 
and UJ 2 ~ 120 ° above. 

The contour for A02 = 30° and the boundary for the libration of 02 are compared to the contour 
for tun = 5 X 10^ yr and the contours for Xu — and 2.0 in Figure 5d, where the positions of the 
fits I-IV are also indicated. We can see from Figure 5d that the fit I with the minimum which 
is unstable (Fig. 2), is in fact far from the stability boundary. Fit II at 62 = 0.14 and ^2 = 130° is 
the fit with the smallest libration amplitudes of both 0i and 02- The evolution of aj, ej, and 9j for 
2000 yr for fit II with sini = 1 is shown in Figure 6 . The semiamplitude A0i ~ 6 ° and A 02 ~ 10°, 
and the small amplitudes ensure that both aj and ej have little variation. Fit III at 62 = 0.02 and 
a ;2 = 130° is an example of a stable fit with large amplitude librations of both 0i and 02 - The 
semiamplitude A0i k. 68 ° and A02 ~ 82°, the eccentricity ei varies between 0.16 and 0.40 and 
62 between 0.02 and 0.33 (see Fig. 7). Finally, fit IV at 62 = 0.02 and a ;2 = 260° is an example 
of a stable fit with only 0i librating. Figure 8 shows that 02 circulates through ±180° when 62 is 
very small. The semiamplitude A0i ~ 94°, the eccentricity ei varies between 0.08 and 0.39 and 62 
between 0.008 and 0.35. 

While the fits with small amplitude librations of both 0i and 02 should be indefinitely stable, 
stability for 5 x 10^ yr does not in general guarantee stability on much longer timescales. Indeed the 
changes in the contours for fun = 2000, 10^, and 5 x 10^ yr in Figure 5a indicate that the region of 
fits with only 0i librating is likely to be significantly smaller for stability on much longer timescales. 
We confirm that the fit IV with only 0i librating, as well as fits II and III, are stable for at least 
10® yr by performing direct integrations of these fits for 10® yr. However, direct integrations of the 
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fits at 62 = 0.02 and W 2 = 360° and 310° show that these fits become unstable after 6.0 x 10® and 
3.4 X 10^ yr, respectively. 

In Figure 9 we present the results of the direct numerical orbit integrations for sinf = 0.5. 
The range of fits that are stable for at least 5 x 10^ yr is significantly smaller than that for sin i = 1 
(compare Figs. 5a and 9o). There are small regions where the fits have only Q\ librating about 0° 
and they are stable for 5 x 10^ yr (see Fig. 9d), but the changes in the contours for tun = 2000, 
10^, and 5 x 10^ yr in Figure 9a indicate that these fits are unlikely to be stable on much longer 
timescales. Thus, for sinf = 0.5, only fits with both 9i and 02 librating about 0° are stable. The 
unstable fit I is again far from the stability boundary. Fit IV, which is stable for at least 10® yr 
for sini = 1, is unstable for sini = 0.5. Direct integrations for 10® yr show that fits II and III with 
sinf = 0.5 are stable for at least 10® yr. We can see from Figures 56 and 96 and Figures 5c and 9c 
that there are only small changes in the contours of the semiamplitudes A 0 i and A 02 for the fits 
that are stable for both sini = 0.5 and 1. This means that the changes in A0i and A 02 with sini 
are small for most of the fits that are stable for both sini = 0.5 and 1. In particular, fit II is also 
the fit with the smallest libration amplitudes of both 9i and 02 for sini = 0.5. Figure 10 (compared 
to Fig. 6 ) shows that the libration period is shorter due to the larger planetary masses but that 
A 01 and A 02 are within 1° of the values for sini = 1. We do not show the time evolution of 9i and 
02 for fit III with sini = 0.5, but A0i and A02 are again within 1° of the values for sini = 1. 


5. DISCUSSION AND CONCLUSIONS 

We have analyzed the HD 82943 planetary system by examining the best-fit double-Keplerian 
model to the radial velocity data as a function of the poorly constrained eccentricity and argument 
of periapse of the outer planet’s orbit. We have not found any good fits that are non-resonant 
and dynamically stable (if the orbits are assumed to be coplanar), and the two planets in the HD 
82943 system are almost certainly in 2:1 mean-motion resonance. The fit I with the minimum 
is unstable, but there is a wide range of fits outside the minimum with only slightly higher than 
the minimum. If the unknown sini ~ 1, there are stable good fits with both of the mean-motion 
resonance variables, 0i = \i — 2 A 2 -|- voi and 02 = \i — 2 A 2 -k ^ 2 , librating about 0° (e.g., fits H 
and HI), as well as stable good fits with only 0i librating about 0° (e.g., fit IV). If sini w 0.5, only 
good fits with both 0i and 02 librating about 0° (e.g., fits H and HI) are stable. Fit H is the fit 
with the smallest libration semiamplitudes of both 0 i and 02 -, with A0i ~ 6° and A 02 ~ 10°. 

Our analysis differs from that of Ferraz-Mello et al. (2005) in having the additional Keck 
data, in using instead of RMS as the primary measure of the goodness of fit, and in examining 
systematically the dynamical properties of all of the fits found. Although Ferraz-Mello et al. (2005) 
showed ej and u)j of many fits with RMS within ~ 0.2 ms“^ of the minimum RMS, they discussed 
in detail the dynamical properties of only two fits — their solutions A and B — with RMS within 
0.06 ms“^ of the minimum. Their solutions A and B have tci ~ 120° and ta 2 ^ 200° and hence large 
amplitude librations of both 0i and 02 - Because the fit to the CORALIE data with the minimum 
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RMS (which is near the minimum Xu of Mayor et ah 2004 at 62 = 0.18 and L 02 = 237°; see Table 
1) is close to the stability boundary in the 62-^2 plane, Ferraz-Mello et al. (2005) were able to find 
stable hts with RMS within 0.06 ms“^ of the minimum. For the combined data set, the fit with 
the minimum RMS at 62 = 0.22 and 002 = 275° is far from the stability boundary, and none of our 
fits with RMS within 0.06 ms“^ of the minimum are stable (see Figs. 36, 5a, and 9a). However, by 
examining systematically the Xu^ RMS, and dynamical properties of all of our fits, we have found 
fits like fit IV (which has stable, large amplitude libration of only 9i if sini ~ 1.0), fit III (which has 
large amplitude librations of both 9i and @ 2 ) and, in particular, fit II (which has small amplitude 
librations of both 9i and 92 ), all with RMS within 0.13ms“^ of the minimum and Xu within 0.08 
of the minimum. The stable and unstable regions in the 62 cos W 2 -e 2 sina ;2 plane shown in Figures 
4 and 8 of Ferraz-Mello et al. (2005) are qualitatively similar to those in our Figures 5a and 9a in 
the 62 - 1^2 plane, but it should be noted that their figures show simulations of their solution B (for 
the CORALIE data) with 62 and uj 2 changed, while our figures show simulations of the fits (to the 
combined data set) that minimize Xu given 62 and uj 2 - 

The relatively large values of Xu 1.84) and RMS (> 7.87 ms“^) of the fits presented in this 
paper suggest that the double-Keplerian model may not fully explain the radial velocity data of HD 
82943. Also hinting at the same possibility is the increase in the RMS of the best fit from 6.99 m s“^ 
(for the extracted CORALIE data alone) to 7.88 ms“^ with the inclusion of the Keck data, which 
fill in some gaps in the CORALIE data and increase the time span of observations. However, the 
RMS and xt values are by themselves not very strong evidence, because our estimate (4.2 ms“^) 
for the radial velocity jitter has an uncertainty of ~ 50% and the jitter can be ~ 6 ms“^. On the 
other hand, there appears to be systematic deviations of the radial velocity data from the double- 
Keplerian fits. Eigure 4 shows that the data within about ±200 days of JD 2452600 are higher than 
the fits and that the data near JD 2453181 are lower than the fits. The fact that ~ 40% of the 
Keck measurements fall in these two regions result in the aforementioned increase in the RMS of 
the best fit from 6.99ms“^ for the CORALIE data alone to 7.88 ms“^ for the combined data set. 
But it is important to note that the Keck and CORALIE data are consistent with each other where 
they overlap and that they are both higher than the fits in the region around JD 2452600. One 
possible explanation for the deviations is that they are due to the mutual gravitational interaction 
of the planets. However, unlike the GJ 876 system where the resonance-induced apsidal precession 
rate dwj/dt is —41°yr“^ on average and the precession has been seen for more than one full 
period (Laughlin et al. 2005), the average apsidal precession rates of, e.g., fits H and HI are only 
—0?51 (sin i)~^ yr~^ and —0?71 (sin i)~^ yr~^, respectively, and the orbits have precessed a negligible 
3-4 (sin i)~^ degrees over the 6.1 yr time span of the available observations. On the other hand, there 
are small variations in the orbital elements on shorter timescales due to the planetary interaction 
that could produce deviations from the double-Keplerian model. Alternatively, the deviations could 
be due to the presence of additional planet(s). Continued observations of HD 82943 combined with 
dynamical fits should allow us to distinguish these possibilities. 
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Table 1. Orbital Parameters of the HD 82943 Planets from Mayor et al. (2004) 


Parameter 

Inner 

Outer 

P (days) 

219.4 

435.1 

K (ms“^) 

61.5 

45.8 

e 

0.38 

0.18 

to (deg) 

124 

237 

M (deg) 

357 

246 


Note. — The parameters are 
the orbital period P, the veloc¬ 
ity amplitude K, the orbital eccen¬ 
tricity e, the argument of periapse 
cn, and the mean anomaly Ai at 
the epoch JD 2451185.1. 
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Table 2. Measured Velocities for HD 82943 from Keck 


JD 

(-2450000) 

Radial Velocity 
(ms“^) 

Uncertainty 

(ms“^) 

2006.913 

43.39 

3.2 

2219.121 

23.58 

2.7 

2236.126 

30.29 

2.7 

2243.130 

38.46 

2.6 

2307.839 

-46.33 

3.2 

2332.983 

-11.26 

3.4 

2333.956 

-12.25 

3.3 

2334.873 

-0.65 

3.0 

2362.972 

25.02 

3.3 

2389.944 

47.81 

3.1 

2445.739 

55.26 

3.0 

2573.147 

-48.70 

2.9 

2575.140 

-46.26 

2.6 

2576.144 

-50.28 

3.1 

2601.066 

-22.60 

3.1 

2602.073 

-15.76 

2.9 

2652.001 

19.60 

3.1 

2988.109 

-89.93 

2.8 

3073.929 

-4.31 

3.1 

3153.754 

0.41 

2.8 

3180.745 

-62.93 

2.7 

3181.742 

-54.04 

3.1 

3397.908 

-99.57 

2.8 
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Table 3. Orbital Parameters of the HD 82943 Planets 


Parameter 

Inner 

Outer 

Fit I Ix^ = 1.87, 

xliw) 

= 1.84, RMS = 7.88ms-i] 

P (days) 

219.3 

441.2 

K (ms“^) 

66.0 

43.6 

e 

0.359 

0.219 

a; (deg) 

127 

284 

M (deg) 

353 

207 

msini (Mj) 

2.01 

1.75 

Vk (ms"^) 


-6.3 

Vc (ms“^) 


8143.9 

Fit II ixl 

= 1.92, 

RMS = 8.00 ms“^) 

P (days) 

219.6 

436.7 

K (ms“^) 

53.5 

42.0 

e 

0.422 

0.14 

uj (deg) 

122 

130 

M (deg) 

357 

352 

msini (Mj) 

1.58 

1.71 

Vk (ms“^) 


-6.9 

Vc (ms“^) 


8142.7 

Fit III ixl 

= 1.91, 

, RMS = 7.99ms-i) 

P (days) 

219.5 

438.8 

K (ms“^) 

58.1 

41.7 

e 

0.398 

0.02 

UJ (deg) 

120 

130 

M (deg) 

357 

356 

msini (Mj) 

1.74 

1.71 

Vk (ms“^) 


-6.3 

Vc (ms“^) 


8143.2 











Table 3—Continued 


Parameter 

Inner 

Outer 

Fit IV {xl 

= 1.90, 

RMS = 7.97ms-i) 

P (days) 

219.5 

439.2 

K (ms“^) 

59.3 

41.7 

e 

0.391 

0.02 

a; (deg) 

121 

260 

M (deg) 

356 

227 

msini (Mj) 

1.78 

1.72 

Vk (ms“^) 


-6.2 

Vc (ms“^) 


8143.3 


Note. — Fits II-IV are 10 parameter fits with fixed 
e and uj of the outer planet’s orbit. Fit I is a 12 pa¬ 
rameter fit, and its Xu shown both at its true value 
and at the value rescaled to 10 adjustable parameters 
for comparison with the 10 parameter fits. In addition 
to the parameters P, K, e, lv, and A4, the minimum 
planetary mass, msinf, and the zero-point velocities, 
Vk and Vc-, of the Keck and CORALIE data, respec¬ 
tively, are listed. 
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Fig. 1.— Radial velocity curve of the double-Keplerian fit I (with the minimum compared to 
the Keck {filled circles) and CORALIE {open triangles) data for HD 82943. The error bars show 
the quadrature sum of the internal uncertainties and estimated stellar jitter (4.2ms“^). 



Fig. 2.— Evolution of the semimajor axes, ai and 02 , and eccentricities, ei and 62 , for fit I with 
sini = 1 and starting epoch JD 2451185.1. The system becomes unstable after ~ 3200 yr. 
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Fig. 3.— Parameters of the best-fit double-Keplerian models for HD 82943 for fixed values of the 
eccentricity, 62 , and argument of periapse, iaJ 2 , of the outer planet’s orbit, (a) Contours of Xu ™ 
steps of 0.05 from 1.85 to 2.0 and then in steps of 0.25 to 3.0. The positions of fits I-IV are marked. 
( 6 ) Contours of the RMS of the residuals (in units of ms“^) in steps of 0.1 from 7.9 to 8.2 and then 
in steps of 0.6 to 10. (c) Contours of the eccentricity, ei, of the inner orbit in steps of 0.02 from 
0.3 to 0.5. (d) Contours of the argument of periapse, wi, of the inner orbit in steps of 10° from 80° 
to 150°. (e) Contours of the planetary mass ratio mijm 2 in steps of 0.1 from 0.7 to 1 . 2 . The ratio 
mi/m 2 is nearly independent of sini, as long as sini is not very small and mi, m 2 <C mo. 



1500 2000 2500 3000 3500 

JD-2450000 

Fig. 4.— Residual velocities from fits I-IV of the Keck {filled circles) and CORALIE {open triangles) 
data for HD 82943. The error bars show the quadrature sum of the internal uncertainties and 
estimated stellar jitter (4.2ms“^). 
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(a) (b) 



Fig. 5.— Dynamical properties of the best-fit double-Keplerian models shown in Fig. 3 for sini = 1. 
(a) Contours of the time tun (in units of yr) at which a model becomes unstable, (b) Contours 
of the semiamplitude A0i in steps of 20°. (c) Contours of the semiamplitude A 62 in steps of 20° 
from 20° to 100°. The thick solid contour is the boundary for the libration of 02 - (d) Positions of 
fits I-IV in the e 2 -cu 2 plane. The gray contours are xt ~ and 2.0. The thick solid contour is 
^un = 5 X 10^ yr. The thin solid contours are A 62 = 30° and the boundary for the libration of 62 - 
Only fits with A 6 i < 140° are stable for at least 5 x 10^ yr. Fit I at the minimum of is unstable. 
Fit II is the stable fit with the smallest libration amplitudes of both 61 and 02- Fit III is a stable 
fit with large amplitude librations of both 61 and 02 - Fit IV is a stable fit with only 0\ librating. 
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Fig. 6 .— Evolution of the semimajor axes, ai and eccentricities, ei and 62 , and 2:1 mean-motion 
resonance variables, 6*1 = Ai — 2 A 2 -|- wi and 62 = \i — 2 A 2 -|- W 2 , for fit II with sini = 1. The 
semiamplitude A 0 i 6 ° and A 02 ~ 10 °. 



Fig. 7.— Same as Fig. 6 , but for fit III with sini = 1. The semiamplitude A0i ~ 68 ° and 
A 6 I 2 ~ 82°. 







circulates. 
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Fig. 10.— Same as Fig. 6, but for fit II with sinf = 0.5. The semiamplitudes A0i and A02 are 
within 1° of the values for sini = 1. 





